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ABSTRACT 

Theories of galaxy formation predict the existence of extended gas halo around spiral 
galaxies. If there are 10-100 nG magnetic fields at several ten kpc distances from the 
galaxies, extended galactic cosmic ray (CR) haloes could also exist. Galactic CRs can 
interact with the tenuous hot halo gas to produce observable gamma-rays. We have 
performed search for a gamma-ray halo around the M31 galaxy - the closest large spiral 
galaxy. Our analysis of almost 7 years of the Fermi LAT data revealed the presence of a 
spatially extended diffuse emission excess around M31. The data can be fitted using the 
simplest morphology of a uniformly bright circle. The best fit gave a 4.7cr significance 
for a 0.9° (12 kpc) halo with a photon flux of ~ (3.2 ± 1.0) x 10 -9 cm _2 s _1 and a 
luminosity of (4.Oil.5) x 10 38 erg s _1 in the energy range 0.3-100 GeV. Our results also 
imply a low level of the flux from the disc of the M31 galaxy (3.3il.O) x 10 _1 ° cnU 2 s _1 . 
The corresponding gamma-ray luminosity, 5 x 10 3 ' erg s _ 1 is several times smaller than 
the luminosity of the Milky Way. This difference could be explained by a lower star 
formation rate in M31: there are less CRs and the level of the ISM turbulence is lower, 
which in turn leads to a shorter time of CR containment. 
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1 INTRODUCTION 


Theories of galaxy formation predict the existence of 
extended haloes around spiral gal axies due to gas in¬ 
flow from their neighb ourghood (IWhite &; Reesl 1 19781 : 
iFukugita fc Peebles! 1200a ). When falling, this gas could be 
heated up to virial temperatures 10 6 — 10 7 K, producing 
huge reservoirs of hot gas (coronae). There are several ob¬ 
servational manifestations of these coronae: soft diffuse X- 
ray emission e xtendin g up to several ten kpc from the 
central galax y (Li! et alJ l2008h, absorption in O VII line 
dWang et al.l 120051 : Bregman fc LIovd-Davie3 SoO^distor- 
tions in the shape of gas clouds ( Westmeier et al . 20(1) and 
stripping of gas in the satellite galaxies by th e ram-pressure 
of th e halo gas dBlitz &: RobishawlkOOOl) . see dPutman et al.l 
and references therein for an extensive review. The 
existence of such a hot halo around th e Milky Way is es¬ 
tablis hed by several different methods dMiller fc BregmarJ 
l2013h . 
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Extragalactic hot haloes are extremely elusive: they 
have been observed only i n several normal ga laxies with¬ 
out star-formation bursts dBogdan et alJ 120131 ). However, 
recently the implementation of stacking technique showed 
that there is a statistically significant excess in the extended 
X-ray signal from som e normal late-type galaxies as well 
dAnderson et al.l I2013h . suggesting that even if the haloes 
are individually undetectable at the present level of sensi¬ 
tivity, they still could be discovered by analysis of groups of 
objects. 

The Milky Way and other disc galaxies can also be 
immersed into extended cosmic ray s (CRs) haloes. This 
idea was thoroughly investigated in dDe Paolis et al.l 1 19991 : 
iFeldmann et al.ll2013l ). It is well known that the Milky Way 
is not a perfect calorimeter for CRs: they rather quickly, 
on time scales of 10-20 Myr, escape from dense regions of 
the Galaxy, losing only minor part of their energy in in¬ 
terac tions with the interstellar medium dStrong et al.ll20071 . 
l20ld) . However, if strong enough magnetic fields exist far 
away from the central regions of the Galaxy, these CRs 
would not go directly away to the intergalactic space, but 
would be instead retained in the magnetized galactic halo 
for a considerable time. Magnetic fields 10-100 times as weak 






































2 M. S. Pshirkov, V.V. Vasiliev & K. A. Postnov 


as the galactic ones 0(fj,G) could be sufficient to contain 
these CRs for the cosmological time. Wandering CRs would 
interact with tenuous (~ 10 -4 cm -3 ) hot plasma produc¬ 
ing gamma-rays via pionic channel. Estimates show that 
the gamma-ray luminosity of such a halo could be around 
10 39 erg s _1 at energies above 100 MeV (iFeldmann et al.l 
1201311 . The size and shape of the halo cannot be firmly es¬ 
tablished and depend crucially on the propagation proper¬ 
ties of CR s. The halo ’half-light ’ radius is estimated to be 
20-40 kpc dFeldmann et alj|2013h . 

The contribution of the CR halo around our Galaxy to 
the isotropic gamma-ray background can be as high as 10%, 
and it is difficult to disentangle it from the truly extragalac- 
tic component. However, such haloes can be searched for 
around other spiral galaxies. The most natural target is the 
halo around the nearby M31 (Andromeda) galaxy. With the 
expected angular size of several degrees and a gamma-ray lu¬ 
minosity of ~ 10 39 erg s _1 , such a halo could be detected by 
the Fermi LAT even from the Earth-M31 distance of > 700 
kpc. The presence of a hot gas around M31, which is essential 
for the gamma-ray emission from the CR halo, was recently 
demonstrated by the discovery of certain absorption features 
in UV-s pectrum of quasars projected on the sky close to the 
galaxy (iRao et al.ll2013l : Lcl iner et a l.l l2014l ) and distortions 
in the observed CMB spectrum i n the vicinity of M31 due 
to interference from the halo gas dDe Paolis et al.||20l3 ). 

The paper is organized as follows: in Section II we de¬ 
scribe the data and method of data analysis, Section III con¬ 
tains our results, and summary and discussion are in Section 
IV. 


2 DATA AND DATA ANALYSIS 

In our analysis we have used 81 months of Fermi LAT data 
collected since 2008 Aug 04 ( MET =239557417 s) until 
2015 Jul 06 (MET=457860004 s). We have selected events 
that belong to the ’’SOURCE” class in order to have a suf¬ 
ficient number of events without losing in their quality. The 
PASS8_V2 reconstruction and vlOrOpqj version of the Fermi 
science tools was used. As the expected signal is weak and 
diffuse, we have selected events with energies larger than 300 
MeV, because at lower energies the Fermi LAT point spread 
function (PSF) quickly deteriorates. Usual event quality cut, 
namely that the zenith angle should be less than 100° (which 
is sufficient at these energies) has been imposed. 

Smaller PSF allowed us to use smaller region of inter¬ 
est (Rol) as well - we took a circle of 10 degrees around 
the centre of the M31 galaxy (aj 2 ooo = 10.6846°, <5J 2000 = 
41.2692°). The data were analysed using the binned maxi¬ 
mum likelihood approach dMattox et al.lll996l ) implemented 
in the gtlike utility, in which two model hypotheses were 
compared by their maximal likelihoods with respect to the 
observed photon distribution. The null hypothesis does not 
include the halo, the alternative hypothesis adds the halo to 
the list of sources of the null hypothesis. 

The source model includes 26 sourc e s from the 3FGL 
catalogue iThe Fermi-LAT Collaboration! d2015l l. the latest 
galactic interstellar emission model giLiem_v06_revl.fit, and 
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the isotropic spectral template iso_source_vO6.tx10. Param¬ 
eters of these sources were allowed to change. We also in¬ 
cluded additional 69 point-like gamma-ray emitters from the 
3FGL catalogue between 10° and 15° from the Rol center 
with their parameters held fixed. 

The M31 galaxy itself was modeled as an 
extended source b ased o n th e IR observations 
dMiville-Deschenes fc Lagachel i2005l l (100/irn normal¬ 

ized IRIS map) fo llowing the presc riptions of the Fermi 
LAT collaboration (lAbdo et alj|201oh . 

Finally, extended halo spatial templates were inserted 
into the source model. We have used the simplest spatial 
models - uniformly bright circles of different radii (from 0.1 
to 5.0 degrees with 0.1 degree step). Of course, it is not 
a realistic model, because some decrease in surface bright¬ 
ness towards the outer halo regions can be expected. On 
the other hand, scarcity of the data used justifies this sim¬ 
ple approach - a more sophisticated model would inevitably 
involve a larger number of parameters, which would make 
fitting much harder and would dilute any obtained signifi¬ 
cance as well. 

The M31 galaxy and the halo were described by a simple 
power-law model: 

dN/dE = N 0 (E/E 0 )~ r (1) 

The normalization No and spectral index P were allowed 
to vary during the likelihood optimisation, while the energy 
scale Eo was fixed at 1 GeV. 

The evidence of the detection of gamma-ray signal from 
the halo was evaluated in terms of a likelihood ratio test 
statistic: 

TS = - 2 In Lmax ’° ( 2 ) 

where Lmax ,0 and Lm aX} 1 are maximum likelihood values 
obtained from the observed data fit using null and alterna¬ 
tive hypothesis, respectively. If the alternative hypothesis is 
true, then y/TS is approximately equivalent to the source 
detection significance. 


3 RESULTS 

Firstly, we have performed our analysis without additional 
source. The M31 galaxy was modelled in two different ways: 
as a point-like source or as an extended object (the IRAS 
template). The extended template for the M31 galaxy fits 
the data considerably better than the simple point-like 
source (TSWt = 79, TS pa = 62.3). The galaxy has a soft 
spectrum with photon index T = 2.40 ± 0.12 and the flux 
F = (2.6 ± 0.4) x 10 _9 ph cm” 2 s _1 in the 0.3-100 GeV en¬ 
ergy range. The spectrum is even softer if the galaxy is mod¬ 
elled as a point-like source: T = 2.64±0.15 with the photon 
flux F = (1.9 ± 0.3) x 10 _9 ph cm -2 s _1 . 

The results of fitting with additional halo component 
are presented in Fig. [T] The fit quality improvement can 
be easily seen. The highest statistical significance TS = 22 
was obtained for a halo with radius Rhaio = 0.9°, corre¬ 
sponding to a linear size of ~ 12 kpc. The photon flux 

2 http://fermi.gsfc.nasa.gov/ssc/data/ 
access/lat/BackgroundModels.html 
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from the extended halo and the 0.3-100 GeV luminosity ob¬ 
tained from the fit are ~ (3.2 ± 1.0) x 10~ 9 cm^ 2 s'* and 
(4.0 ± 1.5) x 10 38 4 erg s _1 , respectively, adopting the dis¬ 
tance d = 780 kpc. The spectral index is found to be rather 
soft: T = 2.30 ± 0.12. A marginal improvement (TS ~ 8 
) could be also achieved by adding a 3-degree halo (~ 35 
kpc). Note that in the case of the small halo the total sig¬ 
nal from the M31 region is dominated by the halo rather 
than the galaxy disc: the flux from the disc is found to be 
~ 10% of the total flux (~ (3.3 ± 1.0) x 10~ 10 cm _2 s _1 ). 
This fact suggests that the IR-based template cannot fully 
trace the gamma-ray emission, and this emission is far more 
extended than the template size. To find how the observed 
gamma-ray flux from the region is shared between the two 
components, we have simulated events in the energy range 
0.3-100 GeV for the relevant time span (71 months) and 
the Rol described above. The model included the galactic 
and isotropic backgrounds, 24 point-like sources from the 
3FGL catalogue, the M31 galaxy disc (taken in the form of 
the IRAS 100/im template). Spectral parameters and pho¬ 
ton fluxes for the point-like sources were taken from the 
3FGL catalogue. The recommended values were taken for 
the isotropic and galactic backgrounds fluxefl We have fixed 
the total flux from the halo and disc components to a fidu¬ 
cial value 3.0 x 10~ 9 cm _2 s _1 and performed simulations, 
gradually changing the halo contribution from 0 to 100 %. 
The results are presented in Fig0 First of all, note that 
there is no leakage of the disc photons to the halo - when 
the fraction of the simulated halo photons is low, the results 
of the corresponding fit immediately show this. The same 
is true for the disc component as well. Clearly, this method 
can effectively separate the two components. 

The actual flux from the M31 galaxy disc is low, < 
4 x 10 -10 cm” 2 s _1 , therefore its luminosity in the 0.3-100 
GeV energy range is rather modest: 5 x 10 3 ' erg s _1 . This 
value is several tim es as small as that of the Milky Way 
(IStrong et al.ll201(il H. This difference could be explained by 
low level of the star formation rate in M31, which is smaller 
than the corresponding Galactic rat e by a factor of 4-5 
dKennic.utt fe Evandhoui ; iFordl l2013l f . The lower star for¬ 
mation rate not only implies less energy in the form of cos¬ 
mic rays that eventually produce the observed gamma-rays, 
but also means that the level of turbulence in the ISM of 
M31 is lower, which in turn leads to a rapid escape of CRs 
from the M31 disc region. 

To exclude possible systematic and instrumental effects, 
which could affect our results, we have performed several 
additional tests: 

(i) In order to check whether the TS increase corre¬ 
sponding to the 0.9-degree halo was caused by some uniden¬ 
tified point-like sources, not included into the 3FGLcata- 
logue, we have calculated the TS map using the gttsmap 
utility (see fig. 0 . A TS exccess at about 0.9 degree from 
the center of the galaxy with the galactic coordinates (l = 
120.58°, b = —21.17°) emerges that could be ascribed to 
FSRQ B3 0045+013. However, even after adding the source 


3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/help/ 
gtobssim.txt 

4 This luminosity also includes some contribution from the Milky 
Way halo, so the direct comparison is not straightforward 



Figure 1 . TS/Tthalo) curves: an earlier version (65 month of 
data, Pass7 Reprocessed events and 2FGL catalogue) is shown 
for comparison. The TS(Rhalo) curve is much smoother when the 
latest version of event reconstruction and 3FGL source catalogue 
are used. 
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Figure 2. Comparison of the reconstructed and simulated fluxes 
from the disc and halo components for different halo fractions in 
the total flux. 

with these coordinates into our source model, the TS of the 
halo decreased only from 22 to 15 (the TS for this source 
was 11.2), thus the whole increase could not be attributed 
to this source alone. Alternatively, this TS excess could be 
produced by an inhomogeneity in the M31 halo. The plau¬ 
sibility of this scenario is also confirmed by the inspection 
of the TS maps of several simulated haloes - they are far 
from being smooth and uniform, but rather consist of several 
random knots that could have TS > 10 (see Fig. 0. 

(ii) We have also checked that the smallness of our Rol 
does not considerably affect our analysis: we have performed 
the data analysis using a larger circle with 15° radius. The 
TS values from the halo remained essentially unchanged. 


4 SUMMARY AND CONCLUSIONS 

Using almost 7 years of the Fermi-LAT data, we have per¬ 
formed searches for an extended gamma-ray halo at ener¬ 
gies above 300 MeV around the closest large spiral galaxy, 
M31. Such a gamma-ray halo could have appeared as a re¬ 
sult of interactions of CRs from the M31 galaxy with gas 
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Figure 3. Templates of the M31 galaxy disc and the halo with 
radius R=1.0°. 



Figure 4. The TS map with the IRAS template for the M31 disc. 
A complex extended structure around the galaxy is clearly seen. 
1- and 3- degree radius white circles are shown for convenience. 



Figure 5. The TS map of simulated data including 0.9-degree 
halo with the photon flux Fo.3— lOOGeV — 1.5 X 10 9 cm 1 . A 
bright spot with TS ~ 10 that emerged by chance is seen. There 
are no clear signs of any extended structure beyond the 1° radius. 


in its halo. We find that the Fermi-LAT data suggest the 
presence of a spatially extended gamma-ray excess around 
M31. The data can be described using the simplest morphol¬ 
ogy of a uniformly bright circle. The best fit gave ~ 4.70- 
significance for a 0.9° radius (12 kpc) halo with the pho¬ 
ton flux ~ (3.2 ± 1.0) x 10 -9 cm _2 s _1 and luminosity 
(4.0 ± 1.5) x 10 38 erg s _1 in the energy range 0.3-100 GeV. 
The presence of such a halo compcllingly shows that a sub¬ 
stantial magnetic field (> 100 nG) should extend around 
M31 up to at least 10 kpc. This excess could indicate the 
presence of a compact CR halo of 10-15 kpc in radius, similar 
to the halo that is indirectly observed around the Milky Way. 
Independent observational checks of such a circumgalactic 
magnetic fields could be done, for example, by analysis of 
the Faraday rotation measures of backgrou nd extragalac- 
tic sources (see, e.g., dPshirkov et al.l 1201 lh 1 . In addition, 
synchrotron emission from secondary leptonic CR compo¬ 
nent could contribute at radio-frequencies ^100 MHz. Fi¬ 
nally, gamma-ray signatures of dark matter particle anni ¬ 
hilations (or decays) aro und M31 may be expected dBaltd 
l2008i : iDugeer et alJl20fo ). 

Our results also imply a low level of the gamma-ray flux 
from the M31 galaxy disc - (3.3 ± 1.0) x 10 _1 ° cm~ 2 s _1 . 
The corresponding gamma-ray luminosity, 5 x 10 37 erg s _1 , 
is several times smaller than the corresponding gamma-ray 
luminosity of the Milky Way. This difference could be ex¬ 
plained by a lower star formation rate in M31: there are less 
CRs and the level of the ISM turbulence is lower, which in 
turn leads to a shorter time of containment. 

Past activity of the M31 galaxy could have been respon¬ 
sible for the complex structure of the TS excess at several 
degrees scale - see the example of the Fermi bubbles in the 
Galactic center. The SMBH in the center of M31 that is al¬ 
most by two orders of magnitude more massive than SMBH 
in the Milky Way and could have injected much more en¬ 
ergy in the form of cosmic rays in the circumgalactic s pace. 
Futu r e observatio ns, including at energies > 100 GeV (IBirdl 
120151 : ISmithll2015l f would certainly clarify this issue. 
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